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C° DISCONTINUOUS GALERKIN FINITE ELEMENT METHODS FOR 
SECOND ORDER LINEAR ELLIPTIC PARTIAL DIFFERENTIAL 
EQUATIONS IN NON-DIVERGENCE FORM* * * § 

XIAOBING FENGt, LAUREN HENNINGS*, AND MICHAEL NEILAN§ 

Abstract. This paper is concerned with finite element approximations of W 2,p strong solutions of second- 
order linear elliptic partial differential equations (PDEs) in non-divergence form with continuous coefficients. A 
nonstandard (primal) finite element method, which uses finite-dimensional subspaces consisting globally contin¬ 
uous piecewise polynomial functions, is proposed and analyzed. The main novelty of the finite element method is 
to introduce an interior penalty term, which penalizes the jump of the flux across the interior element edges/faces, 
to augment a nonsymmetric piecewise defined and PDE—induced bilinear form. Existence, uniqueness and error 
estimate in a discrete W 2,p energy norm are proved for the proposed finite element method. This is achieved by 
establishing a discrete Calderon-Zygmund-type estimate and mimicking strong solution PDE techniques at the 
discrete level. Numerical experiments are provided to test the performance of proposed finite element method 
and to validate the convergence theory. 


1. Introduction. In this paper we consider finite element approximations of the following 
linear elliptic PDE in non-divergence form: 

(1.1a) Cu := —A : D 2 u = f in 

(1.1b) u = 0 on dfl. 


Here, D C 


is an open bounded domain with boundary <9D, / G L p (yt) (1 < p < oo) is 


given, and A = A(x) £ [C°(n) 


is a positive definite matrix on £7, but not necessarily 


differentiable. Problems such as (1.1) arise in fully nonlinear elliptic Hamilton-Jacobi-Bellman 


equations, a fundamental problem in the field of stochastic optimal control mm- In addition, 
elliptic PDEs in non-divergence form appear in the linearization and numerical methods of fully 
nonlinear second order PDEs 

Since A is not smooth, the PDE (1.1a) cannot be written in divergence form, and therefore 
notions of weak solutions defined by variational principles are not applicable. Instead, the 
existence and uniqueness of solutions are generally sought in the classical or strong sense. In 
the former case, Schauder theory states the existence of a unique solution u £ C' 2,a (D) to 
(1.1) provided the coefficient matrix and source function are Holder continuous, and if the 
boundary satisfies <9£7 £ C 2 ’ a . In the latter case, the Calderon-Zygmund theory states the 
existence and uniqueness of u £ IE 2,P (£7) satisfying (1.1) almost everywhere provided / £ L p (£7), 
A £ [C'°(0)] rixn and <9£7 £ C 1 ’ 1 . In addition, the existence of a strong solution to (1.1) in two- 
dimensions and on convex domains is proved in mm®. 

Due to their non-divergence structure, designing convergent numerical methods, in par¬ 
ticular, Galerkin-type methods, for problem (1.1) has been proven to be difficult. Very few 
such results are known in the literature. Nevertheless, even problem (1.1) does not naturally 
fit within the standard Galerkin framework, several finite element methods have been recently 
proposed. In m the authors considered mixed finite element methods using Lagrange finite 
element spaces for problem Of An analogous discontinuous Galerkin (DG) method was pro¬ 
posed in (9j. The convergence analysis of these methods for non-smooth A remains open. A 
least-squares-type discontinuous Galerkin method for problem (1.1) with coefficients satisfying 
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Fig. 1.1. Outline of the convergence proof. 


the Cordes condition was proposed and analyzed in |23j . Here, the authors established optimal 
order estimates in h with respect to a H 2 - type norm. 

The primary goal of this paper is to develop a structurally simple and computationally 
easy finite element method for problem (1.1). Our method is a primal method using Lagrange 
finite element spaces. The method is well defined for all polynomials degree greater than one 
and can be easily implemented on current finite element software. Moreover, our finite element 
method resembles interior penalty discontinuous Galerkin (DG) methods in its formulation and 
its bilinear form, which contains an interior penalty term penalizing the jumps of the fluxes 
across the element edges/faces. Hence, it is a C° DG finite element method. In addition, we 
prove that the proposed method is stable and converges with optimal order in a discrete W 2,p - 
type norm on quasi-uniform meshes provided that the polynomial degree of the finite element 
space is greater than or equal to two. 

While the formulation and implementation of the finite element method is relatively simple, 
the convergence analysis is quite involved, and it requires several nonstandard arguments and 
techniques. The overall strategy in the convergence analysis is to mimic, at the discrete level, 
the stability analysis of strong solutions of PDEs in non-divergence form (see [HI Section 9.5]). 
Namely, we exploit the fact that locally, the finite element discretization is a perturbation of 
a discrete elliptic operator in divergence form with constant coefficients; see Lemma 3.1 The 


first step of the stability argument is to establish a discrete Calderon-Zygmund-type estimate 
for the Lagrange finite element discretization of the elliptic operator in <o> with constant 
coefficients, which is equivalent to a global inf-sup condition for the discrete operator. The 
second step is to prove a local version of the global estimate and inf-sup condition. With these 
results in hand, local stability estimate for the proposed C° DG discretization of (1.1) can be 
easily obtained. We then glue these local stability estimates to obtain a global Garding-type 
inequality. Finally, to circumvent the lack of a (discrete) maximum principle which is often used 
in the PDE analysis, we use a nonstandard duality argument to obtain a global inf-sup condition 
for the proposed C° DG discretization for problem (1.1). See Figure for an outline of the 
convergence proof. Since the method is linear and consistent, the stability estimate naturally 
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leads to the well-posedness of the method and the energy norm error estimate. 

The organization of the paper is as follows. In Section [2] the notation is set, and some 
preliminary results are given. Discrete W 2,p stability properties, including a discrete Calderon- 
Zygmund-type estimate, of finite element discretizations of PDEs with constant coefficients are 
established. In Section[3j we present the motivation and the formulation of our C° discontinuous 
finite element method for problem (1.1). Mimicking the PDE analysis from [14] at the discrete 
level, we prove a discrete W 2 - p stability estimate for the discretization operator. In addition, 
we derive an optimal order error estimate in a discrete W 2 ’ p -norm. Finally, in Section |4j we 
give several numerical experiments which test the performance of the proposed C° DG finite 
element method and validate the convergence theory. 


2. Notation and preliminary results. 

2.1. Mesh and space notation. Let II C M" be an bounded open domain. We shall 
use D to denote a generic subdomain of D and dD denotes its boundary. W S ’ P (D) denotes the 
standard Sobolev spaces for s > 0 and 1 < p < oo, W°- P (D) = L P (D) and Wg’ p (n) to denote 
the subspace of W S,P (Q) consisting functions whose traces vanish up to order s — 1 on dSl. (v)d 
denotes the standard inner product on L 2 (D) and (•, •) := (-, -)q. To avoid the proliferation of 
constants, we shall use the notation a < b to represent the relation a < Cb for some constant 
C > 0 independent of mesh size h. 

Let Th '■= 7/t(H) be a quasi-uniform, simplical, and conforming triangulation of the domain 
Q. Denote by £■{ the set of interior edges in Th , the set of boundary edges in Th, and 
£h = £{ U , the set of all edges in Th- We define the jump and average of a vector function 
w on an interior edge e = dT + D dT~ as follows: 


H 
M 

where w ± = ia| T± and n± is the outward unit normal of T ± . 

For a normed linear space A', we denote by X* its dual and (•, •) the pairing between X* 
and A'. The Lagrange finite element space with respect to the triangulation is given by 

(2.1) V h := {v h € Hq(Q.) : v h \ T € P fc (T) VT G %}, 

where P*,(T) denotes the set of polynomials with total degree not exceeding k (> 1) on T. We 
also define the piecewise Sobolev space with respect to the mesh Th 

W s,p (Th) := I] WS ’ P ( T )> W h ] ■= W 2 ' p (Th) H Wq 1 ’^), 
tgu 

L P h {T h ) := n L p (T ), W°’ P (D) := W s ’ p (T h )\ D , L P (D) := L p (T h )\ D - 

TeT h 


= ■ n_J 


w 


• n_ 


1 ( , 

= - l W 1 • 71 + — W 

2 V +| - 


n_ 




For a given subdomain D C H, we also define 14(D) C 14 and W^\d) C W ( h P> as the subspaces 
that vanish outside of D by 

V h (D) : = {ve V h - u|n\ c = 0}, W (p \d) := {v G W {p) - v\ n \ D = 0}. 

We note that 14(D) is non-trivial for diam(D) > 2 h. 

Associated with D C SI, we define a semi-norm on W^' P (D) for 1 < p < oo 

(2.2) M\ w ^ d) = II D 2 v\\ lhd) + ( Y, ^IIMH^^nS)) 1 - 
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Here, D^v £ L 2 (Cl) denotes the piecewise Hessian matrix of v, i.e., D 2 v\t = D 2 v\t for all 
T£%. 

Let Qh : L p (Cl) —> Vh be the L 2 projection defined by 
(2.3) (Qhw,v h ) = (w,v h ) \/w G L 2 (Cl), v h £ V h . 

It is well known that [8] Qh satisfies for any w £ W m ’ p (Cl) 


(2.4) || Qhw\\w^p(n) ^ IMI w^p(n) m = 0,1; 1 < p < oo. 

For any domain D C Cl and any w £ L p h {D), we also introduce the following mesh-dependent 
semi-norm 


(2.5) 




Ll{D) ■= . SUp 


(w,V h ) 


D 


0 ?v h EV h (D) \\ v h\\LP'(D) 


where —I-= 1. 

P P' 


By (2.3), it is easy to see that || • || L p h (D) is a norm on Vh{D). Moreover by (2.4) 

( 2 . 6 ) 


i M ( Wh,v) ( w h ,Q h v 

\wh\\LP(n) = sup TT - TT1 -= sup 


< 


vglp'( n) IPllLp^n) veLp'(n) ll w lli>'(n) 

(w h , Q h v) 


sup 

eLp'(n) II 2^ l, llLp'(n) 


< H^IU^n) Vwh £ Vh- 


(v) 

2.2. Some basic properties of Wf functions. In this subsection we cite or prove 

(p) 

some basic properties of the broken Sobolev functions in W h , and in particular, for piecewise 
polynomial functions. These results, which have independent interest in themselves, will be 
used repeatedly in the later sections. We begin with citing a familiar trace inequality followed 
by proving an inverse inequality. 

Lemma 2.1 ([!]). For any T £ Th, there holds 

(2-7) |M| P L p {dT) < (^ _1 ||Vw||^ (T) + h^\\v\\l P{T) ) V« e W 1,P (T) 

for any p £ (l,oo). Therefore by scaling, there holds 


( 2 . 8 ) 


E 

ed£ 


IIP _ < 

iLP(eri-D) ~ 


h 


I INI Um 
IIHIL (0) 


h p \\Vv\\ p LP{D) 


to £ V h {D), 
\/v £ W^ P \D). 


Lemma 2.2. For any Vh £ Vh, D C Cl, and 1 < p < oo, there holds 


(2.9) 

\\ v h\\wf’ p (D) ~ h 1 \\ v h\\w 1 -p(D h ) 

where 


(2.10) 

D h = {x £ Cl : dist(x, D) < h} 

Proof. By (2.2), (2.7) and inverse estimates |7J|I], we have 


\\ v h\\wl' p (D) ~ \\^h v h\\LP(D) + ( E ? ll ILp(enD)) 

eesl 
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<\\D 2 h v h \\ LP{D )+ E (4- p (^r 1 pviii P( T ) +^ 1 iiv^r iP(T) )) ,; 

TGTfe 

TCD h 


□ 

The next lemma states a very simple fact about the discrete W 2 ’ p norm on W^’ p (fl). 
Lemma 2.3. For any 1 < p < oo, there holds 

(2.11) IMI wf ,p (n) — IMlwMfi) € W ,P (H). 


Next, we state some super-approximation results of the nodal interpolant with respect to 
the discrete W 2,p semi-norm. The derivation of the following results is standard m , but for 
completeness we give the proof in Appendix [A] 

Lemma 2.4. Denote by Ih : C°(f2) —> Vh the nodal interpolant onto Vh- Let rj £ C°°(D) 
with fo d~ J for 0 < j < k. Then for each T £ Th with h < d < 1, there holds 

(2.12) h m \\pv h - Ih(vv h )\\ Wm , P{D ) < ^\\vh\\Lp(D h ) for m = 0,1, 

(2-13) \\r)v h - Ih(vv h )\\ w 2 , P{D) < -p\\v h \\ w i,r( Dh ), 

Moreover, if k > 2, there holds 

(2-14) || Vv h - h{riv h )\\ w 2 , P{D) < ^\\v h \\ W 2 , P{Dh) . 

Here, D C Dh C satisfy the conditions in Lemma \2.2\ 

To conclude this subsection, we state and prove a discrete Sobolev interpolation estimate. 
Lemma 2.5. There holds for all 1 < p < oo, 

H v HlL(n) ^ IMUp(n)IMIwt p ( fi) e w h ] - 


|V^f iP( o) 


(2.15) 


Proof Writing ||Vw||^ p(n) = f n |Vw| p 2 Vw ■ S7wdx and integrating by parts, we find 

J (jVu>| p ~ 2 Afc + (p— 2)\\7w\ p ~ 4 (D 2 w\7w) ■ '\7w S jwdx 
+ E / [[|Vw| p " 2 Vw]]u;ds 

E J \Vw\ p ~ 2 \Df l w\\w\ dx + E J [[|Vw| p ~ 2 Vu;]]u;ds. 


T&T h 


eesl 


To bound the first term in (2.15) we apply Holder’s inequality to obtain 

(2.16) ^jVH p - 2 |£ 2 HMdz< HlV^I^I^^HD^IUp^lkllL^o) 
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Likewise, by Lemma |2.1| we have 


(2.17) 


J [[|Vw| p 2 Vtv^wds 

E (5 II MIL ne)T~ 2 { h ^ llMlU’(e))(^|MU*>(e)) 
Z l|Vw|li; ( n)lkllip(n)( E ^MlMDIlEe))" 


e££ 

< 


eesl 


~ ll v HIL^)IHURn)IMIw^(fi)- 

Combining (2.15 )-(|2.17|) we obtain the desired result. The proof is complete. □ 


2.3. Stability estimates for auxiliary PDEs with constant coefficients. In this 


subsection, we consider a special case of (1.1a) when the coefficient matrix is a constant matrix, 
A{x) = A 0 £ R nx ”. We introduce the finite element approximation (or projection) £ 0 ,/t of A on 

(v) 

Vh and extend A,ft to the broken Sobolev space W h . We then establish some stability results 
for the operator £o,fc,- These stability results will play an important role in our convergence 
analysis of the proposed C° DG finite element method in Section [3j 
Let A q £ R nx " be a positive definite matrix and set 

(2.18) Cqw := — A 0 : D 2 w = —V • (A q \7w). 

The operator C 0 induces the following bilinear form: 


(2.19) ao(w, v ) := (Cqw, v) = / AoVw-Vvdx Vw, v £ Hq (fl), 

Jn 

and the Lax-Milgram Theorem (cf. jTD]) implies that £q 1 : H~ 1 ( 17) — > 7Lg(f2) exists and is 
bounded. Moreover, if dfl £ C 1,1 , the Calderon-Zygmund theory (cf. [T3J Chapter 9]) infers 
that £q 1 : L p (fit) —> W 2 ’ p (fl) fl W 0 1,p (fl) exists and there holds 

(2-20) HA) Vllw 2 .p(fi) < IMlLp(fi) Vtp £ L p (fl). 

Equivalently, 

(2-21) IMk 2 -(n) < || A^|Lp(o) € W2 ’ P n 

The bilinear form naturally leads to a finite element approximation (or projection) of A 
on Vh, that is, we define the operator £ 0 ,fc, : V), —> V),, by 

(2.22) (C 0 ,hw h ,v h ) := a 0 {w h ,v h ) \/v h ,w h £ V h . 


Remark 2.1. When A = I, the identity matrix, £o,h is exactly the finite element the 
discrete Laplacian that is, A,ft = — A^. By finite element theory [3], we know that A, ft : A —► 
Vh is one-to-one and onto, and therefore : Vh —> Vh exists. 

Recall the following DG integration by parts formula: 


(2.23) 


Vhvdx 


/ (y h ■ t) v dx + 
Jn 


Y (/MM 

eGfl e 


ds 
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/ {{ T }} ' H ds ) + /( T ' n e)vds, 

Je ee£? Je 


which holds for any piecewise scalar-valued function v and vector-valued function r. Here, 


is defined piecewise, i.e., V/Jt = V |t for all T £ Th■ For any Wh,Vh £ 14, using (2.23) with 
r = Ao'Vwh , we obtain 


(2.24) a 0 (w h ,v h ) = - / (H 0 : Dlw h )v h dx + ^ 

Jn 


v h ds. 


We note that the above new form of ao(-, ■) is not well defined on i?o(H) x Uq(H). However, 
it is well defined on Wjk x wjf 1 with f + T = 1- Hence, we can easily extend the domain 
of the operator £ 0 ,/i to broken Sobolev space Wff\ Precisely, (abusing the notation) we define 

£o,h ■ —> (Wff ^)* to be the operator induced by the bilinear form ao(-.-) on W^ x W ^ \ 

namely, 

(2.25) (jC 0 ,hW,v) := a 0 (w,v) Vw £ , v £ W^ p } . 

A key ingredient in the convergence analysis of our finite element methods for PDEs in 
non-divergence form is to establish global and local discrete Calderon-Zygmund-type estimates 


similar to (2.21) for These results are presented in the following two lemmas. 

Lemma 2.6. There exists h o > 0 such that for all h £ (0, ho) there holds 


(2.26) 


I w h || wf ,p (Q) ~ Vwft £ 14- 


Proof. First note that (2.26) is equivalent to 
(2.27) 


£ 


o,hTh\\w^ p (n) ~ llWilk(n) £ 14- 

For any fixed ip h £ 14, let w := £ IU 2 ’ P (H) n W 0 1,P (H) and w h := £ k- 

Therefore, w and Wh, respectively, are the solutions of the following two problems: 


(2.28) 


a 0 {w,v) = (tp h ,v) Vu£iLg(H), a 0 {w h ,v h ) = (<p h ,v h ) \/v h eV h , 


and thus, Wh is the elliptic projection of w. 
By (2.21) we have 


(2.29) 


Mlrv 2 .p(n) ~ ll^llip(n)- 


Using well-known L p finite element estimate results [4j Theorem 8.5.3], finite element interpo¬ 
lation theory, and (2.29) we obtain that there exists ho > 0 such that for all h £ (0, ho) 


(2.30) \\w — uih\\w 1 ’i’(n) ^5 k — -4Hlwi'P(n) < ^kk 2 ' p (fi) ^5 h\\tph\\LP(n)- 

It follows from the triangle inequality, an inverse inequality (see Lemma|2.2[), the stability 


of Ih , (2.29) and (2.30) that 


Ik - w h\\w^ p ( fi) ~ k “ J hW\\wl' p ( n) + 114^ - w h\\wf p (Q) 

^ k’k 2 .p(n) + h^Whw - w h || ff i,p( fl) 
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Thus, 


||<P/i||i>(fi) + h 1 \\IhW — w\\ W i, P (n' ) + h 1 ||u> — Wft|| l yi,p(f 2 ) 


\\ w h\\ w ?^<si) ~ ll w - w h\\ w ^(n) + IMIvkMo) ~ ll^ll^(n), 


which yields (2.271, and hence, (2.261. □ 

Lemma 2.7. For xq G 12 and R > 0, define 

(2.31) Br(x o) := {x £ 12 : \x — £o| < i?} C 12. 

Let R' = R + d with d > 2h. Then there holds 

(2-32) \\ Wh \\wl' p (B r (x 0 )) ~ 11^0, hWh\\L^(B R , (x 0 )) ^ w h £ Vf,(Bji(x o)), 


Proof. To ease notation, set Br := Br(x o) and Ur, := -Br'(:to). Recalling (2.6), we have 
by Lemma |2~G| 

\\ W h\\w?’ p (B R ) = \\ W h\\w^(n) ~ l|A,^|| LP (a) < ||£o,/ l U'fc||i*(n) = sup 

h K ’ h y ’ v h ev h F/i|lip'(n) 

Set R" = (R + R')/2, so that R < R" < R 1 . Denote by Xb r „ the indicator function of 
Bru := Br"{xq). Since % = 0 on LI\Br, we have 

a 0 {w h ,v h ) = a 0 (w h ,XB R „v h ) = a 0 (w h , Ih(XB R „v h )) Vv h £ V h . 

Moreover, we have Ih(XB R „v h ) € V h (B R >) and 

\\Ih(XB Rl ,Vh)\\Lp' (b r ,) = \\Ih(XB Rl ,Vh)\\Lp' (Q) < IIXBr// v h || lp’ (f2) < \\ v h Wlp' (fj) • 
Consequently, 


a 0 {w h ,I h (xB R „v h )) a 0 {w h ,v h ) 

\\ w h\\w?’ p (B R ) - SU P TiTT-yi| -< sup tt— n - 

h ’ v h ev h \\Ih(XB R „v h )\\ LP ’ {B ) Vh ev h (B Rl ) \\ v h\\ L p'( b rI ) 


= sup 


(£ 0 ,hWh,Vh) 


v h £V h (B R ,) \\ v h\\Lp’(B R/ ) 


- = \\£o,hWh\\L p h (B Rl )- 


□ 


3. C° DG finite element methods and convergence analysis. 

3.1. The PDE problem. To make the presentation clear, we state the precise assump¬ 
tions on the non-divergence form PDE problem B Let A £ [C'°(D)]" X " be a positive definite 
matrix-valued function with 


(3.1) 


A|£| 2 A(:r)£-£< A|£| 2 V(eR B ,ieH 


and constants 0 < A < A < oo. Under the above assumption, C is known to be uniformly elliptic, 
hence, strong solutions (i.e., IU 2 ’ P solutions) of problem (1.1) must satisfy the Aleksandrov 
maximum principle mmm- 










By the W 2,p theory for the second order non-divergence form uniformly elliptic PDEs m 
Chapter 9], we know that if dtt £ C 1,1 , for any / £ L P (Q) with 1 < p < oo, there exists a 
unique strong solution u £ W 2,p (fl) fl W^’ P (D.) to satisfying 


(3.2) 


Mlw^.jqn) ^ ||/||z,p(f 2 )- 


Moreover, when n = 2 and p = 2. it is also known that ns m o Em m the above conclusion 
holds if fl is a convex domain. 

For the remainder of the paper, we shall always assume that A £ [C 0 (fi)] raxn is positive 


Calderon-Zygmund estimate (3.2). 


definite satisfying (3.1), and problem (1.1) has a unique strong solution u which satisfies the 


3.2. Formulation of C° DG finite element methods. The formulation of our C° DG 

finite element method for non-divergence form PDEs is relatively simple, which is inspired by the 
finite element method for divergence form PDEs and relied only on an unorthodox integration 
by parts. 

To motivate its derivation, we first look at how one would construct standard finite element 
methods for problem (1.1) when the coefficient matrix A belongs to [C ,1 (fl)] nxra . In this case, 


since the divergence of A (taken row-wise) is well defined, we can rewrite the PDE (1.1a) in 
divergence form as follows: 


(3.3) 


-V • (AWu) + (V • A) • Vu = /. 


Hence, the original non-divergence form PDE is converted into a “diffusion-convection equation” 
with the “diffusion coefficient” A and the “convection coefficient” V • A. 

A standard finite element method for problem (3.3) is readily defined as seeking Uh £ 14 
such that 

(3.4) [ (AVu h ) -Vv h dx + f (V • A) -Vu h v h dx = j fv h dx Mv h £V h . 

JCl J Q 

Now come back to the case where A only belongs to [C' 0 (f2)]" xn . In our setting, the 
formulation (3.4) is not viable any more because V- A does not exist as a function. To circumvent 
this issue, we apply the DG integration by parts formula ( 2.23[ ) to the first term on the left-hand 
side of (3.4) with r = AViih and V in (3.41 is understood piecewise, we get 


(3.5) - / (A : Dlu h )v h dx + ^2 iA'Vu h §v h ds= / fv h dx Vv h £V h . 

Jn eg£j r Je J n 

Here we have used the fact that [W]] = 0 and Vh\an = 0. 

No derivative is taken on A in (3.5), so each of the terms is well defined on 14. This indeed 


yields the C° DG formulation of this paper. 

Definition 3.1. The C° discontinuous Galerkin (DG) finite element method is defined by 
seeking Uh £ 14 such that 


(3.6) 

where 


a h {u hl v h ) = (f,v h ) Vv h £V h , 


(3.7) a h (w h ,v h ):=- (A : P 2 h w h )v h dx + ^ 

Jn 


v h ds, 
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(3.8) 


( f,Vh) ■= [ fv h dx Vv h £V h . 
J n 


A few remarks are given below about the proposed C° DG finite element method. 
Remark 3.1. (a) The above method is also defined for A G [L°°(fl)] I, ' xrl and no a priori 
knowledge of the location of the singularities of A are required in the meshing procedure. 


(b) The C° DG finite element method (3.6) is a primal method with the single unknown 
Uh■ It can be implemented on current finite element software supporting element boundary 
integration. 


(c) From its derivation we see that (3.6) is equivalent to the standard finite element method 


(3.4) provided A is smooth. In addition, if A is constant then (3.6) reduces to 

a 0 (u h ,v h ) = (f,v h ) \/v h £V h . 

This feature will be crucially used in the convergence analysis later. 

(d) In the one-dimensional and piecewise linear case (i.e., n = 1 and k = 1), the method 


(3.6) on a uniform mesh equivalent to 

A(xi)( - Ci -1 + 2 a - c i+ i) = h 2 f(xi), 

where Uh = Y2iLi > an d \}P'h}f=i represents the nodal basis for Vh- 

3.3. Stability analysis and well-posedness theorem. As in Section |2.3[ using the 
bilinear form a^(-, •) we can define the finite element approximation (or projection) C /, of C on 
Vh, that is, we define £ h : Vh —> Vh by 


(3.9) 


(C h w h ,v h ) := a h (w h ,v h ) \/v h ,w h G V h . 


Trivially, (3.6) can be rewritten as: Find Uh G V) ( such that 

(£ h Uh,v h ) = (f,v h ) Vv h G V h . 


Similar to the argument for £q_ h, we can extend the domain of Ch to the broken Sobolev space 
W^\ that is, (abusing the notation) we define Ch : W^ 


(Wl p,) )* by 


(3.10) 


(C h w,v) := a h (w,v) Vw G W^\ v G W } 


Ap) 


The main objective of this subsection is to establish a W 2 ' p stability estimate for the 
operator Ch on the finite element space V),. From this result, the existence, uniqueness and 


error estimate for (3.6) will naturally follow. The stability proof relies on several technical 


estimates which we derive below. Essentially, the underlying strategy, known as a perturbation 
argument in the PDE literature, is to treat the operator Ch locally as a perturbation of a stable 
operator with constant coefficients. The following lemma quantifies this statement. 

Lemma 3.1. For any S > 0, there exists Rs > 0 and hg > 0 such that for any Xo G Cl with 
A) = A(x$) 


(3.11) 


\{C h - C 0 ,h)w\\ L P h{BRsixo)) < <5|Mlw^ ( b k .(xo)) Vw g W h ] ' Wh ^ hs - 


Proof. Since A is continuous on fl, it is uniformly continuous. Therefore for every 5 > 0 
there exists R$ > 0 such that if x,y G D. satisfy \x — y\ < Rs , there holds |A(x) — A{y)\ < 6. 
Consequently for any x 0 G SI 

m - Ao|| L oo (B j?i5 ) < 5 
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(3.12) 









with B Rs := B Rs (x 0 ). 

Set hs = min{/io, and consider h < hs, w £ Wjf' 1 and Vh £ Vh(B Rs ). Since (£o,fc — 
A)™ £ W^ p) , it follows from (Of, ([2724]), (O), ([322]), and Of that 


((A),h - Ch)w,v h ) 

= - [ ((A) - A : D 2 h w)v h dx+ ^ [ 1(A) - A Vio] v h ds 

JBr s eg£ / AnB^ 

< ||2- A||L~(B J?i5 )(||£ ) ^|| L p( BRj )||u h || Lp /( Bji ^ 

+ (£ ft £ 2 1l IV,,.]] r„ ( , ns „,)' (£ A) 

ee££ ee£^ 

< ||2 — AIU~(b Bi 5 )IAII^, P ( B )||wh|| iP '( B ) < ^ll w ll^’ p (s B )II w?i IIlp'(_b B s )- 


The desired inequality now follows from the definition of || • \\l^(b r )■ D 
Lemma 3.2. There exists R\ > 0 and hi > 0 srtc/i that for any xq £ O 


(3.13) 


II^AI wA( b ri( X o)) 1$ IA/i' w ftllLg(B Ji2 (Eo)) S Vh(Bf j 1 (a;o))) V/i ^ A> 


with R 2 = 2Ri. 


Proof. For d 0 > 0 to be determined below, let R\ = \Rs 0 as in Lemma 3.1 Let hi = % 


and set Bi = B Ri (x o). Then by Lemmas 2.7 and 3.1 with d = Ri and A = A^o); we have for 
any w h £ V h (Bi) 

ll W ^llrV^ p (B 1 ) ~ \\£o,hW h \\ L P {B2 ) < ||(£o,/» - £h)Wh\\ L P(B 2 ) + \\£hWh\\L°(B 2 ) 

1$ ^olAh.ll W^ ,F (B 2 ) 4“ I \£-hWh\\L p h (B 2 ) = ^olA^II w^’ p (Bi) + W^-hWh || L p (B 2 ) • 

For <5 0 sufficiently small (depending only on A), we can kick back the first term on the right-hand 
side. This completes the proof. □ 


Lemma 3.3. Let Ri and hi be as in Lemma 3.2. For any xq £ Cl, there holds 
(3.14) ||AH| L £(b R i(xo )) ~ IHIw^ p (B Rl( x o)) \/w ewl p) ,\/h<hi. 


Proof Set Bi = B Rl {x o). By the definition of Ch, (2.6), (2.8) and (2.4), we have for any 
Vh £ V h (Bi) 


{C h w,v) = - (A : Dlw)v h dx + ^2 / [[AVw]]u h ds 

J B\ £~c i " gX~\B\ 

etc. 


< 


IAaIU^boIMI 

B p, ( B i) 


< I \D 


(E h l 1l[[ VM, ]]llB P (enB 1 )) P (E M«hll £,'( enSl 

ee£l ee£‘ h 

(B 1 ) + (E^- P ||M||L(e)) P )lKII^(B 1 ) 


>h w \\L* 


e6 £1 


\wf ,p (Bi) 11 Vh 11 Lp 1 (Bj) • 
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The desired inequality now follows from the definition of || • || B P( Bl ). □ 

Lemma 3.4. Let hy be as in Lemma \3.°2\ Then there holds for h < hy 

(3-15) \\ w h\\ w ^p(Q) < \\£hWh\\LP(n) + ||^||iP(n| Vwft £ V h . 


Proof. We divide the proof into two steps. 

Step 1: For any Xq G Q, let Ry and hy be as in Lemma [372] let i ?2 = 2Ry, R 3 = 3i?i, and 
set By = B Bi ( 20 ) for i = 0,1,2. Let 77 £ C 3 (fl) be a cut-off function satisfying 


(3.16) 


0< V <1, r?| Bi = l, V\ n \ B2 =0, 


We first note that rjWh £ wI p \b 2 ) and Ihfnwh) £ 14 (B 3 ) for any Wh £ 14. Therefore, by 
Lemmas 2.4 (with d = Ry) and 3.2 we have 

\\ W h\\wf' p (B 1 ) = \\ 1 T w h\\w£' p (B 1 ) - 1 V w h — Ih(ll w h)\\]y2,P( Bl ) + ||-4 (f] w h) || 


\w—m = 0(d~ m ) m = 0 , 1 , 2 . 


< 


-^2\\w h \\w 1 ’P(B 2 ) + \\Ih(riWh)\\ w 2.r( Bl ) 


~ ^2 \\ W h\\w^P(B 2 ) + \\^h(Ih(l7Wh))\\Ll(B 2 ) 


< 


^2 IW| W i,p(B 2 ) + \\£h(VWh)hl(B 2 ) + II Ch(r]Wh - Ih(VWh))\\Ll(B 2 )- 


Applying Lemmas 3.3 and 2.4 we obtain 

1 


(3.17) 


\\ W h\\wl’ p { Bl ) ~ fi2\\ w h\\w'-P(B 2 ) + II^OwOlU^Bs) 

+ \\ryw h - I h (vw h )\\ w 2 ,P( B2) 

< ^2 ll u; ^llw i -p(B 3 ) + \\c h (vw h )\\ LliB3) . 

To derive an upper bound of the last term in (3.17), we write for Vh £ 14(7?3), 

(C h {rru} h ),v h ) = - l A : D\{ryw h )v h dx + ^ / \[A\7(r)w h )]]v h ds 

J B 3 pczf 1 ^ en ^3 

= -/ (r)A : D\wh + 2 AX7r] ■ Vw h + w h A : Dlrfjvh dx + ^ / [[AVto/,]]^ ds 

cp/ ^ e n^ 3 

= (C h w h , I h (r]v h )) - / (2AV?? • Vwh + : D 2 h rf)v h dx 

Jb 3 

- {A: Dlw h )(r]Vh-Ih(vvh))dx+ ^2 / [[AVtu/J (^ _ hivvh)) ds. 

J B 3 ct-ei J eDB 3 


By Holder’s inequality, Lemmas 2.1 2.2 2.4| and (3.16) we obtain 

3—21 


{C h {r]W h ),V h ) < ||£7 l MJ/ l |Up (B 3 ) ||4(m)ILp' (B3 ) +^ 1 2 |kfc|kl.J>( B 3)ll^lli I >'( B 3) 


+ \\ w h\\w^ p (B 3 ) 


(ll VVh - Ih(vvh)\\ LP ’ { B 3 ) + h\\V(vvh - IhirjVh)) || t p'(B 3 )) 
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< (\\£ h Wh\\ L F h ( B3 ) + ^2 II w ^I|vF1>p(B3)) \\ v h\\LP'(B 3 'j^ 

which implies that 

\\£h{riWh)\\Ll(B 3 ) ^ \\£hW h \\L p h (B 3 ) + ~j^\\wh\\w^P(B 3 )- 


Applying this upper bound to (3.171 yields 
(3.18) 


ll^llw^C-Bi) ~ \\£hWh\\L*(B 3 ) + ^2 \\ w h\\w l ’P(B 3 ) VlCfc, € V h . 


Step 2: We now use a covering argument to obtain the global estimate (3.15). To this 
end, let C f l with N = 0(R^[ n ) sufficiently large (but independent of h ) such that 

f l = UjL 1 Bn 1 (xj). Setting Sj = Br 1 (xj) and Sj = Bn 2 (xj) = B 2 R 1 {xj), we have by (3.18) 


N 




3 = 1 


P < 

wZ' p (Sj 


N 


3=1 

N 


~ J^2p ll^hllvy 1 '®^) + ll^^lli 


3 = 1 


p „ + — 
L U S >) Rl p 

LlCSjY 




Since Vh{Sj) C Vh, we have 


N N 

E 11^^11^^) = E 

t=i 


3= 

< N 


(C h w h ,v h ) 
sup i|„ II 

II 

(c h w h ,v h ) 
sup l„ II 

0?v h eV h (Sj) WVhWLp'iSj) 

3=1 

0 ^v h £V h (Sj) IrClLp'(fi) 


(C h w h ,v h ) 


sup 

o ^tvh&Vh ll u ft|lLp'(n) 


< 


R r { 


\\£hWh\\ P L P(ny 


Consequently, since i?i is independent of h, we have 


\ w h\\w?' p (Q.) ~ -^\\£hWh\\L p ((i) + -52 ll u ’^lliv 1 .!>(n) ^ \\£hWh\\iP{n) + ||wi?i||iv 1 .p(n)- 

h p ill 


Finally, an application of Lemma |2.5| yields 

1 1 

\\ w h\\w£’ p (n) ~ Ili^w/illLjKfi) + ll u; ^ll£p(n)ll u '^ll^2,p^ n y 

Applying the Cauchy-Schwarz inequality to the last term completes the proof. □ 

Using arguments analogous to those in Lemma |3.4[ we also have the following stability 
estimate for the formal adjoint operator. Due to its length and technical nature, we give the 
proof in the appendix. 

Lemma 3.5. There exists an h 2 > 0 such that 


(3.19) 


ll^llip'(n) ~ SU P 


(■ £ h w h ,v h ) 


Yv h e V h 


O^Wh&Vh \\ w h\\w^’ p (p,) 


provided h < h* := min{hi, h 2 } and k > 2. 
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Remark 3.2. Denote by C* h the formal adjoint operator of Ch.- 
equivalent to the stability estimate 


Then inequality (3.19) is 


(3.20) 


IKIIlp'(Q) $ sup 

Ojtw h ev h 


C C* h v hl w h ) 

W' w h\\\y%' p (Q) 


Vv h € V h . 


Thus, the adjoint operator C* h is injective on Vh- Since Vh is finite dimensional, C* h on Vh is an 
isomorphism. This implies that Ch is also an isomorphism on Vh ; the stability of the operator 
is addressed in the next theorem, the main result of this section. 

Theorem 3.1. Suppose that h < min{/ii, / 12 }, and k >2. Then there holds the following 
stability estimate: 


(3-21) ll^llw^n) ~ \\£h' w h\\L v h (n) h € V h - 

Consequently, there exists a unique solution to ( |3.6| ) satisfying 
(3-22) ll' u h.|lvE^' p (n) ~ ll/IU p (0)- 


Proof. For a given Wh £ V/,, Lemma [3.5| guarantees the existence of a unique iph £ Vh 
satisfying 


(3.23) 


{ChV h ,ifh) = / w h \w h \ p 2 v h dx \/v h €V h . 
Jn 


By (3.19) we have 


( C h vh,if h ) J n w h \w h \ p2 v h dx ||P _i 

Lp'(n) sup tt— rr -= sup - 77—77 -< \\wh\\LP (nr 

O^VhtVh \\ v h\\wf L ’ p (Q.) 0 jiv h e.V h ll^fcll W^’ P (Q) 

The last inequality is an easy consequence of Holder’s inequality, Lemma |2.5| and the Poincare- 
Friedrichs inequality. Taking Vh = uih in ( 3.23| ), we have 

IMIiP(n) = {^hW h , 1 p h ) < \\C h W h \\ L r {n) \\lph\\ L p' (Cl) < 11^^11^(0)11^11^), 

and therefore 


IIw^HlpCD) ^ \\ChWh\\L p h (Q.). 


Applying this estimate in (3.4) proves (3.21). 

Finally, to show existence and uniqueness of the finite element method (3.6) it suffices to 
show the estimate (3.22). This immediately follows from (3.21) and Holder’s inequality: 




w% ,p (n) ~ ll-£/i u /tll-L£(n) — _ sup 


(C h u h ,v h ) 


o^v h £V h ll v ft|lLp'(n) 


= sup 


In frh dx 


o^tv h ev h ItMIlp^o) 


< ll/llip(fj)- 


□ 
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3.4. Convergence analysis. The stability estimate in Theorem [3TT] immediately gives us 
the following error estimate in the Wff p semi-norm. 

Theorem 3.2. Assume that the hypotheses of Theorem\3.l\ are satisfied. Let u £ W 2,p (Ll) 


and Uh £ 14 denote the solution to (1.11 and (3.6), respectively. Then there holds 
(3.24) ll w — u h\\w?' p (Q) ~ \\ u — w h\\w?' p (Cl)- 

Consequently, if u £ bF s,p (fi), for some s > 2, there holds 

ll M — u h\\wf’ p (n) ~ ^ ll M lliv f ’P(n)) 

where t = min{s, k + 1}. 

Proof. By Theorem |3.1| and the consistency of the method, we have Vvh £ 14 


II u h ~ w h \ 


W, 


2 ’ p (n) ~ 


\\£h(Uh - W h )\\ L *(Cl) = sup 
o^v h ev h 


= sup 

O^Vh&Vh 


= sup 

O^Vh&Vh 


a h {u h -w h ,Vh) 

ll^llzy (n) 
{C h {u - w h ),Vh) 
I II ip' (0) 


= sup 

0 ^Vh&Vh 

< lb 


( Ch(u h - w h ),v h ) 

ll^llLp^n) 

a h {u - w h ,v h ) 


\ v h\\Lp'(Q) 


u - w h\ 


Wf p (Q)- 


Applying the triangle inequality yields (3.24). □ 


4. Numerical experiments. In this section we present several numerical experiments to 
show the efficacy of the finite element method, as well as to validate the convergence theory. 
In addition, we perform numerical experiments where the coefficient matrix is not continuous 


and/or degenerate. While these situations violate some of the assumptions given in Section 3.1 
the tests show that the finite element method is effective for these cases as well. 

Test 1: Holder continuous coefficients and smooth solution. In this test we take 
12 = (—0.5, 0.5) 2 , the coefficient matrix to be 


A = 


rc| x / 2 + 1 
-M 1 / 2 


-4I 1 / 2 

5|ar| x / 2 + 1 


and choose / such that u = sin(27ra;i) sin(7ra:2) exp(xi cos(x 2 )) as the exact solution. 

The resulting H 1 and piecewise H 2 errors for various values of polynomial degree k and 
discretization parameter hare depicted in Figure |4T] The figure clearly indicates that the errors 
have the following behavior: 


W~ Uh\m(Q.) = 0{h L ), \\Df(u — Uh)\\L 2 (£i) = 0{h k x ). 


The second estimate is in agreement with Theorem |3.2| In addition, the numerical experiments 
suggest that (i) the method converges with optimal order in the Ft 1 -norm and (ii) the method 
is convergent in the piecewise linear case (k = 1). 

Test 2: Uniformly continuous coefficients and W 2,p solution. For the second set 
of numerical experiments, we take the domain to be the square LI = (0,1/2) 2 , and take the 
coefficient matrix to be 

5 


/ 


A = 


iog(M) 

l 


15 


15 


1 


1 


iog(M) 


+ 3 , 














Fig. 4.1. The H 1 (left) and piecewise H 2 (right) errors for Test Problem 1 with polynomial degree k = 
1,2, 3,4. The figures show that the H 1 error converges with order 0(h k ), where as the piecewise H 2 error 
converges with order 0{h k ~ 1 ). 


We choose the data such that the exact solution is given by u = |x|'/ 4 . We note that u £ 
W m ’ p {yi) for (7 — 4 m)p > —8. In particular, u £ W 2,p ( fi) for p < 8 and u £ W /3,p (f2) for 
p < 8/5. 

In order to apply Theorem |3.2| to this test problem, we recall that the fcth degree nodal 
interpolant of u with k > 2 satisfies 


h 


for p < 2. Since u £ W 3,p (fl) for p < 8/5, Theorem 3.2 then predicts the convergence 
\\Dl(n-u h )\\ l 2 ( n) < C\\D 2 h (u~I h u)\\ LHn) =0{h 3 /^) 


rate 


for any e > 0. Note that a slight modification of these arguments also shows that |u — 
l h u\ m{n) =0(h 7/4 - £ ). 

The errors of the finite element method for Test 2 using piecewise linear, quadratic and cubic 
polynomials are depicted in Figure 4.2 As predicted by the theory, the H 2 error converges with 
order a; 0{h 3 / 4 ) if the polynomial degree is greater than or equal to two. Similar to the first 
test problem, the numerical experiments also show that the H 1 error converges with optimal 
order, i.e., |u - u h \ H i {n) = 0(/i 7 / 4 “ e ). 

Test 3: Degenerate coefficients and W 2,p solution. For the third and final set of test 
problems, we take If = (0, l) 2 , 


A 


16 ( xf 3 -xf\T 

Q 1/3 1/3 2/3 

9 \ — a; l x 2 X 2 


and exact solution u = xf 3 — xf 3 . We remark that the choice of the matrix and solution is 
motivated by Aronson’s example for the infinity-Laplace equation. In particular, the function u 
satisfies the quasi-linear PDE = 0, where := ( D 2 uS7u) ■ S7u = ( D 2 u ) : (Vtt(Vw) T ). 
Noting that A = Vrt(Vu) T , we see that —A : D 2 u = 0 =: /. 

Unlike the first two test problems, the matrix is not uniformly elliptic, as det(A(cc)) = 0 
for all x £ S2. Therefore the theory given in the previous sections does not apply. We also note 
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Fig. 4.2. The It 1 (left) and piecewise H 2 (right) errors for Test Problem 2 with polynomial degree k = 
1,2,3. The figures show that the If 1 error converges with order fc,7 /4 —where as the piecewise H 2 

error converges with order 0(/i mm 'f fc ’ 7 / 4 — £ 1 ~ 4 ). 


that the exact solution satisfies the regularity u G W m ’ p (Q) for (4 — 3 m)p > — 1, and therefore 
u G W 2 ’ p (ft) n W 1 ’ 00 ^) for p < 3/2. 

The resulting errors of the finite element method using piecewise linear and quadratic 
polynomials are plotted in Figure |4~3l In addition, we plot the computed solution and error in 


Figure 4.4 with k = 2 and h = 1/256. While this problem is outside the scope of the theory, 
the experiments show that the method converges, and the following rates are observed: 


I u - u h || L2(n) = 0(/r 4 / 3 ), |« - u h 1 * 1 ( 0 ) = 0(/i 5/6 ). 



Fig. 4.3. The H 1 (left) and piecewise H 2 (right) errors for the degenerate Test Problem 3 with polynomial 
degree k = 1 and k = 2. The figures show that the L 2 error converges with order ~ (9(/i 4 / 3 ) and the H 1 error 
converges with order ~ 0(/i 5 / 6 ). 
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Computed Solution Error 



Fig. 4.4. Computed solution (left) and error (right) of test problem 3 with k = 2 and h = 1/256. 
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Appendix A. Super approximation result. Here, we provide the proof of Lemma |2.4| 
As a first step, we use standard interpolation estimates ram to obtain for 0 < m < k + 1, 

(A.l) h mp \\rjv h - I h ( V v h )\\ p wm , P(T) < h^\ V v h \ p wk+1 , p{T) . 

Since \r)\wi,°°(T) 1$ d~ 3 and \ v h\H k + 1 (T) = 0, we find 


(A.2) 


\WVh\w k + l *(T) ~ E 

|a:| + |/3|=fc+l ' 
k 


\D a ri\ p \D f> v h \ p dx 


< 


E 




< 


E 


hr™ 


dp(k+l-j) 1^1 Wi’P(T) ~ (p>(k+l-j) \\ Vh \\LP(T)i 

3=0 j =0 


where an inverse estimate was applied to derive the last inequality. Combining (A.21 with (A.l I 
and using the hypothesis h < d then gives us 

*L hP(k+l~j) UP 

h P \\r]Vh — Ih{V v h)\\w™,p(T) ~ E Jp(k+l-i) \\ v h\\ P L*(T) ~ ^ \\ V h lliP(T)‘ 

3=0 


Therefore for m £ {0,1} we have 


h mp \\r]v h - I h (riVh)\\w™,p( D ) ^ E - h{r)v h ) ||^m, P(T) 

T&T h 

TnD^0 


up up 

~ E ^H^ll LP(T) - -^W V h\\LP(D h y 

T&T h 

TnD^d 


Thus, (2.12) is satisfied. 


To obtain the second estimate (2.131, we first use (A.l), (A.2) an an inverse estimate to get 


/j p(fc-i) 


(A.3) \\D~{rjVh Ih^T/Vh)) 11 (T’) ~ h P ^ ' ^ \V v h\y^ k + 1 ,p(T) ~ (p>(k+l—j) I Wi’ p (T) 


\v h \\ 


~ fj2p \\ V h\\ P LP(T) + E dp(k+l-j) IKII W l ’P(T) 

3 = 1 

~ ^2^ll w /»llz,»(T) + ^p\\ v h\\w-i-,p( T ) ~ ( p^ll t; h|l VVUp(t)- 


3=0 

k h k ~ j 


By similar arguments we find 

(A.4) h- p \\\/(r!v h - I h (vv h m P L P {T ) Z ^ k ~%v h f w 


W k+i , P ( T ) ( p.p\\ Vh Ww 1 ’P{T)' 
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Therefore by Lemma 2.7 and dot ( |A.4[ ), we obtain 
\\r]Vh - Ih{vv h )\\ P w 2 , P{D) < Y T h{r]Vh))\\ P LP{T) 


TeT h 

thd^<D 


Y h e P || l^ivvh - Ihivvh))} | 


V 

LP(eCD) 


ee si 


< 


Y [||-D 2 (^-4(m))IIIp(T) + /l - I h (r]Vh))\\ P LP{T) 

T&Th 

TnDyt® 

~ Y, ^ll^llw^-jyT) — ^II^HvKi.p^)- 

TeTh 

TnD^0 


Taking the pth root of this last expression yields the estimate (2.131. The proof of (2.14) uses 
the exact same arguments and is therefore omitted. 


Appendix B. Proof of Lemma 3.5 To prove Lemma 3.5 we introduce the discrete 
fF _2 ’ p -type norm 


(B-l) 


_ _ (r,v h ) D 

I wr 2 ’ p (D) SU P II n 

0^v h eV h (D) \\ v h\\w 2 ’P' (D) 


and the W 1,p norm (defined for L p functions) 

( r,v) D 


(B.2) 


Irll w-^p(d) — 


sup 

O^vZW 1 *' (D) ll u ll W' 1 "*’' (J3) 


sup (r, v) d dx. 

v eff'-r 1 
II IV 1 . p' (£>) — 1 


The desired estimate (3.19) is then equivalent to 
(B.3) 


IKIIzy(fi) ~ \\£'h v h\\ w -2,p' (Qj 'dvh £ Vhi 

where we recall that C* h is the adjoint operator of £/,. Due to its length, we break up the proof 
of (B.3) into three steps. 

Step 1: A local estimate. 

The first step in the derivation of (13.19 ) (equivalently, (B.3)) is to prove a local version of 
this estimate, analogous to Lemma 3.2 To this end, for fixed xo £ D, let <5o, Rs 0 , Ri ■= ^Rs 0 
and B 1 := Bji 1 (x o) be as in Lemmas 3.1 3.2 with <5 0 > 0 to be determined. For a fixed 
Vh £ Vh(Bi), let ip € VF 2,p (fl) D Wq’ p (Q,) satisfy Cip = Vh\vh\ p ~ 2 in f1 with 


(B.4) IMk^n) < lll^r'- 1 |Up ( q) < HW^y 

Multiplying the PDE by Vh, integrating over Q, and using the consistency of Ch yields 


up'- 1 


'Lp’(B i) II h "LP'(Q.) 

Therefore, for any tpn £ Vh.-, there holds 
(B.5) II- 'l p ' 


V h\\ P Tv f fB ^ = IKIIrp'/oN = (A^h) = {C h ip,v h ). 


\ v h\\ P LP ' {Bl ) = (d: h ip h ,v h ) + (C h (<p - iph),v h ) 
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= {C*h v h, Vh) + (£-0,h(V - <Ph), Vh) + ({Ch- Co ,h)(v - <Ph)>Vh), 


where Co,h is given by (2.22) with A 0 = A(x o). Now take iph £ Vh to be the elliptic projection 
of Lp with respect to Co,hi ke., 

(£o,h{<P~ <Ph),Wh) = 0 Vw h £V h . 

Lemma |2.6| ensures that g>h is well-defined and satisfies the estimate 

(B.6) \Wh\\ w ^P(Q.) < \\Co,hVh\\LHp.) = l|A),/i<P||l,*(fi) ^ IMItV 2 .*>(n) ^ \\ V h\\ P Lv '\ Bl y 


Combining Lemma 3.1 (B.4 |-( |B.6[ ) and ( |B.1[ ), we have 

\\ Vh W P Lp'(B,) = ( C h v hiVh) + (( c h - £ 0 ,h)(<P - <Ph),Vh) 

— W~ 2,p ' (fi) ll^^ll Wjj’ p ({2) I \(Ch ~ Co,h)(v — fh)\\L^ i (B 1 )\\ v h\\Lp' (B!) 

— \\Ch Vh \\wp 2 ’ pl (Bi^b^Lp' (B!) — Vh\\ W^' p (Bi) 11^ lllj p/ (Bi) 

< + <5o||wh||^ (Sir 

Taking <5o sufficiently small and rearranging terms gives the local stability estimate for finite 
element functions with compact support: 


(B.7) 


Vh 


\lp'(B 1 ) 1$ \\Ch v h\\ w -2,p' ( Bl ) ^ V h £ Vh(B\). 


Step 2: A global Garding-type inequality. 

We now follow the proof of Lemma |T4| to derive a global Garding-type inequality for the adjoint 
problem. Let i?i be given in the first step of the proof, i? 2 = 2 R ll and R 3 = 3i?i. Let y £ C 3 (fl) 
satisfy the conditions in Lemma [3N| (cf. ( |3.16 )). By the triangle inequality and (B.7) we have 
for any v h £ V h 

IKIlzy (Bi) = llmlliP'(Bi) < Wvvh - hivvh) W L p\ B i ) + \\Ih{r}Vh)\\ L p'{ Bl ) 

< || vv h - Ih(vvh)\\ LP ’ {Bl ) + \\CUIh(vVh))\\ w -2, P ' {Bi) 

< || f]Vh - h(vvh) \\ L p'( B i ) + II C* h {I h (yv h ) - VVh) H^- 2 , p' (B i) + \\C* h (vv h )\\ w - 2 ,p' (Bi) . 


Applying Lemmas 3.3 Lemma 2.4 (with d = R±) and an inverse estimate yields 
( B - 8 ) IKILp'(B i) ~ II Wh ~ Ih{VVh)\\ L p'(B 2) + \\C* h (vvh)\\ w - 2 , P '( Bi) 


< 


< 


^II^IIlp'(b 3 ) + W c h(wh)\\ w -2, P ' {Bs) 


1 


Jji ll^/illw-i'P^Bs) + \\Ch(jl v h)\\ w -2,p' ^ B3 y 


The goal now is to replace C* h {r]Vh) appearing in the right-hand side of (B.8) by C* h v h plus 
low-order terms. To this end, we write for Wh £ Vh(B 3 ) (cf. 2.22), 

(B.9) (. C* h (yv h ),w h ) = a h (wh, yv h ) = a h {w h y, v h ) + [a h (w h , yv h ) - a h (w h r], v h )] 

= a h (I h (w h r]), v h ) + a h (w h y - I h (w h y), v h ) + [a h {w h , yv h ) - a h (w h y, v h )] 

= '■ h + I'2 + lo- 
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To derive an upper bound of I\, we use (B.l) and properties of the interpolant and cut-off 
function rj: 


(B.10) 


= (£* h V h ,I h ( V W h )) < \\£*h v h \\ W - 2 ,p'( B3 )\\I h ( r l w h)\\w 2 P(B 3 ) 


< \\£-h v h\\w- 2 ’V'(B 3 )^ r l w h\\wl’ v (B 3 ) ~ ^2 \\^*h v h || w~ 2 ’P' (B 2 )\\ Wh II W^' P (B 3 ) • 
Next, we apply Lemmas |3.3[ [2A| and an inverse estimate to bound I 2 ' 


(B.ll) I 2 = {C h (r)w h - Ih(ww h )),v h ) < \\r)w h - Ih(vwh)\\ w ^p {B3) \\vh\\ L p'( b 3 ) 

^3 \\ w h\\w^ p (B 3 )\\ V h\\Lp'(B 3 ) < ^3 ll w ftllw^ p (B 3 )ll' l 'ft I!w-Up'(_B 3 )- 


< 


To estimate I 3 , we add and subtract ao(wh , r]Vh ) — oo(w/ l ??, Vh) and expand terms to obtain 

(B.12) / 3 = a 0 {w hl r]v h ) - a 0 (w h r ], v h ) 

+ [ 1 a h (w h ,r]v h ) - a h (w h .i'i, Vh) - (a 0 (w h ,r]v h ) - a 0 (w h ri, v h ))] 

= -J (w h A 0 : D 2 r] + 2A 0 \7ii-\7iUh)vhdx 

- J (w h {A - A 0 ) : D 2 r) + 2(A - A 0 )\/r] • \7wh)v h dx =: Ad + A' 2 . 
Applying Holder’s inequality and Lemmas |C . 1 1 |D. 1 1 yields 

(B.13) Ad < Hu^Ao : D 2 r)\\ w i, P ( B3 )\\v h \\ w -i, P > ^ + 2 / (A 0 X7r] ■ \7w h )v h dx 

J Bs 

~ (^3 IWIwup^) + ^2 ll^llw*- p (B 3 ))IKIIw-Up'(B 3 ) 


< 


\\ w h\\ W £P(B 3 )\\ V h\\w-'’V' (B 3 )-> 


Similarly, by Lemma C.l and (3.12), we obtain 


(B.14) 


K 2 < ||A - A 0 || z ,oo( B3 )(||w; ft || i p( B3 )||D 2 77|| i oo (n ) + || Vu; h || B p( B3 ) || VtjII^^)) ||w h || iP /( S3 ) 
< S 0 (Rl\\D 2 i]\\ L a O ^ n) + i? 3 ||V»7||L~(n))||t«/i|| w -2,P( B3 )||u/»||i,p'( B3 ) 

fio\\ w h\\w£’ p (B 3 ) IMILp' (S 3 ) 


Combining (B.12)—(B.14 ) results in the following upper bound of I 3 : 

(B.15) ^3 1$ ll u ’ /i ll w^ ,p (b 3 ) IK'^II w- 1 ’? 1 (b 3 ) A <5q ||Hd ,p (B 3 ) \\vh || B p'( Ba ) 


Applying the estimates to (B.10) (B.ll), (B.15) to (B.9) results in 


(■ £*h(vv h ),w h ) < -^3 


I W~ 2 'P' (B 3 ) + W'VhWw-i’P 1 (B 3 )' S )\\ w h\\wZ- p (B 3 ) 

^o||||L p'(B 3 ) \\ w h\\ w£' p (B 3 ) ' 
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and therefore by (B.l), 


(B.16) \\Cfr(r]Vh)\\ w -2,p' < ^3 (\\£h V h\\w- 2 ’P' (B 3 ) + II II (. 63 ) • 


Combining (B.16) and (B. 8 ) yields 




1 


I w ~ 2 ’ p '( b 3 ) + II^IIw-i.p'(b 3 ) ) + ^I|%ILp'(b 3 ). 

Finally, we use the exact same covering argument in the proof of Lemma |3.4| to obtain 
II^ILp'(a) ~ W^*h v h\\ w -^,p' ll^llw-cp^a) + ^oII^/iIIlp'^)- 

Taking So sufficiently small and kicking back the last term then yields the Garding-type estimate 

(B-17) ll^frII lp' (st) ~ W^*h v h\\ w -^y ^ + IWIw-^'in)- 

Step 3: A duality argument 


In the last step of the proof, we shall combine a duality argument and (B.17) to obtain the 
desired result (B.3). 

Define the set 


X = {9 S Wq P (Q,) : ||ff||w 1 p(n) = 1}. 

Since X is precompact in L P (Q), and due to the elliptic regularity estimate ||</?||1$ 
||^||LP(fi), the set 

W = {ip £ W 2 ’ p n W 0 1 ,p (f2) : £<p = g, 3g £ X} 

is precompact in W 2,P (Q). Therefore by J22J Lemma 5], for every e > 0, there exists a 
/i 2 (e, W) > 0 such that for each ip £ W and h < hi there exists <ph £ 14 satisfying 


rp 2 .PfO\||'Cfi||z ; P , (Q)- 


1*5.18,) IIV 5 - VhWw^iQ) - £ IOr K - A 

Note that ( |B.18[ ) implies \\Wh\\ w *,p (p.) < ||<^||w 2 .p(n) + £ < 1- 

For g £ X we shall use ip g £ W to denote the solution to Ctp g = g. We then have by 
Lemma [3731 for any vh £ 14 and tph £ 14, 

/ v h gdx = ( c h ip g ,v h ) = (. C* h v h ,ip h ) + (C h {<Pg - Vh),Vh) 

■Jn 

1$ \\£h v h\\ w -*,p' (n)ll^ ,l llw 2 - p (n) W'Aa ~ 7 , ^llw^ ,p (n)l 

Choosing iph so that (B.18) is satisfied (with <p = <p g ) and using the definition of the W~ 1,p 
norm (B.2) results in 

IIDiHh/- 1 ,?' (O) 1$ \\£h V h\\ w -2,p’ (Q) + e ll' i; /i|lLp'(n)- 
Finally we apply this last estimate in ( B.17| ) to obtain 

Ik/illip^a) < \\£'*h v h\\ w -2,p’(p) + e ll^lli,p , (n)' 
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Taking e sufficiently small and kicking back a term to the left-hand side yields (B.3). 
completes the proof. 


This 


Appendix C. A discrete Poincare estimate. 

Lemma C.l. There holds for any Wh £ Vh(D) with diam(D) > h, 


\\wh\\w^p(D) < diam(D) 2 m \\w h \\ w 2, P{D) m= 1,2. 


Proof. Denote by V Ct h C H 2 {fl) n Hq(LI) the Argyris finite element space [4|, and let 
Eh ■ Vh Vc,h be the enriching operator constructed in [5] by averaging. The arguments in [5] 
and scaling show that, for Wh £ Vh(D), 


(C.l) 


E h Wh £ Hy(Dh), 


where Dh is given by (2.101. 
inequality gives 


|| w h EhWh\\ W m ’P(D) fS ll' w h.|lw^ ,p (D) — 0,1,2), 

Since EhWh £ H^Dh) and diam(D) > h, the usual Poincare 


\\E h w h \\ Wm , v{Dh) <diam(D ?i ) 2 m \\E h w h \\ W 2, P{Dh) < diam(D) 2 m \\w h \\ w 2. P{D) . 
Therefore by adding and subtracting terms, we obtain for m = 0,1, 


\\u!h\\w m ’P(D) < \\EhWh\\w m ’P(D) + ~ E h Wh\\w m ’P(D) 

< diam(D) 2 - m \\w h \\ w 2 , P{D) + h 2 - m \\w h \\ w 2 , P(D) 

< diam.{D) 2 - m \\w h \\ w 2 . P{Dy 

where again, we have used the assumption h < diam(D). The proof is complete. □ 

Appendix D. A discrete Holder inequality. 

Lemma D.l. For any smooth function rj, and Wh £ Vh(D),Vh £ Vh, there holds 


(Vry • ViUftK dx 


< 


ID 


I W' 2 ’°°(D)\\ , Wh\\ w z.'P( D )\\Vh\\w~ i y (D)' 


Proof. Let Eh : Vh —t V c be the enriching operator in Lemma C.l satisfying (C.l). Since 
EhWh £ H 2 (D) we have 


id 


(V?y • V (E h w h ))v h dx < || V ?7 • V{E h w h )\\ w i, P ( D ' ) \\vh\\ w -i,pi( D ) 

< \\r]\\w 2 ’°°(D)\\EhWh\\w 2 ’P(D)\\Vh\\ W -i, P ' ( D ) 
W 2 '°° (D) \\ w h\\wf- v (D) W V h\\w- 1 ’P\D) 


< 


Combining this estimate with the triangle inequality, (C.l I, and an inverse estimate gives 

/ (V 77 • S 7 w h )v h dx = / (V;y • V(E h w h ))v h dx + / (V 77 • - E h Wh))v h dx 

J D J D J D 

< \\ri\\w 2 ’°°{D)(\\Wh\\ W 2 ’P( D )\\Vh\\w-l,v'( D ) + II Wh - E h W h \\ W i.P{D)\\vh\\ LP '(D)) 

W 2 >°° (D) \\ w h\\wl’ p (D) ll^hllw- 1 **'^!?)' 


< 
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